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Abstract 

The Ising chain consisting of the antiferromagneticaly coupled ferromagnetic trimer is 
considered in the external magnetic field. In the framework of the transfer- matrix formal- 
ism the thermodynamics of the system is described. The magnetization per site (m) is 
obtained as the explicit function of the external magnetic field (H). The corresponding 
plots of m(H) are drawn. Two qualitatively different regions of the values of coupling 
constants are established: weak antiferromagnetic coupling (Ja < 3Jf) and the strong 
antiferromagnetic coupling (J a > 3Jf). For the latter case the magnetization curve with 
plateau at m/m sa t = 1/3 is obtained. It is proven that the plateau is caused by the stabil- 
ity of spatially modulated spin structure (3111). The values of magnetic field determining 
the width of the plateau are obtained in the limit of zero temperature. 
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1 Introduction 



Quite recently, the family of known non trivial quantum effects in the condensed matter physics 
is enriched with the novel phenomenon — intermediate plateaus in the magnetization processes. 
Namely, the system ceases to responde on changes of the external magnetic field within some 
range of the latter. As a consequence, for some intermediate values of H (H\ < H < H 2 ) a 
horizontal region forms in the magnetization (m) versus the external magnetic field (H) curve 

These plateau structures were experimentally observed in series of compounds: the triangu- 
lar antiferromagnets (CqEu [|l|, CsCuCIj, @, RbFe{MoO<±)2 ||) ; quasi one-dimensional Ni- 
compounds[||]; quasi two-dimensional compound with structure of Shastry-Sutherland lattice 
(SrCu2{BO-z)2 HI). At present there are numerous theoretical papers where the magnetization 



plateaus were obtained for various non-homogeneous Heisenberg chains [)B|-||10||, for spin-ladders 



||1 1|| , and as well as for some two dimensional models [|i2| -||ll 

In the majority of plateau mechanisms which have been proposed up to now the purely 
quantum phenomena play a crucial role. The concepts of magnetic quasiparticles and the 
strong quantum fluctuations are regarded to be of first importance for understanding of these 
processes. Particularly, for a number of systems it was shown that the plateaus at m ^ are 
caused by the presence of the spin gap in the spectrum of magnetic excitations in the external 
magnetic field. Another mechanism lies in so called crystallization of the magnetic particles 



17 



Nevertheless, we suppose that the quasiclassical picture in which the main dynamics of 
spins can be presented in terms of spin "up" and spin "down" is also relevant to the problem 
of magnetization plateaus. The problem formulated in the quasiclassical language allows to 
use the Ising type variable instead of the quantum spin operators. The advantages of this 
approximation are obvious: in the case of one dimensional systems one can develop purely 
analytic and quite simple technique for the magnetization plateau problem. 

One of the pioneering theoretical works where the appearance of the magnetization plateau 
was predicted is the paper of Hida || . In this paper the ferromagnetic-ferromagnetic-anfiferromagnetic 
Heisenberg chain was considered as a model of 3C11CI2 ■ 2dioxane compound which consists of 
the antiferromagneticaly coupled ferromagnetic trimers. The Hamiltonian, considered by Hida 
has the following form: 



H = Y.(^i nm + Ki int ) + Hz, (1) 

i 

where 



Hi trim = _2J F ( Si -T,, + T i -<T i ), 

Hi int = 2J A (cr t - Sl+1 ), (2) 
U z = -gfi B Hj2(s t z + cTi Z +n z ), 

i 

here J a and Jp are the antiferomagnetic and ferromagnetic coupling constants respectively, 
Si, Ti and <Tj are the S = 1/2 spin operators and Hz is the standard Zeeman Hamiltonian. Using 
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the method of the numerical diagonalization of finite size system Hida obtained a magnetization 
curve with the plateau at m/m sat = 1/3, where m sat is the saturation magnetization. 

The approach we are offering in the present paper is based on very simple idea to use 
the Ising spins instead of the Heisenberg operators in Eq. (|l|) and to develop the transfer- 
matrix formalism (see for example Ref. Jl9| ). It is worth to emphasize that in contrast to 
the majority of existing approaches to the problem of magnetization plateau, our technique 
is entirely based on analytical calculations and allows to obtain the magnetization profiles for 
arbitrary finite temperatures and arbitrary values of parameters J a and Jp. As a result we 
obtained the diagrams of magnetization processes for emerged ferromagnetic-ferromagnetic- 
antiferromagnetic Ising chain with the same magnetization plateau as in Ref. ||. 

Presently, there are the numerical results indicating the existence of magnetization plateau 



in two dimensional Ising models [plj. Moreover, in the Ref. |2(J the multisite interaction Ising 
model on the infinite dimensional recurrent lattice was considered. There the magnetization 
curves with plateaus at m = and m/m sat = 1/2 were also obtained there. However, till 
now, there are no any evidence of the appearance of the magnetization plateaus in the one 
dimensional Ising chains. 

2 The ferromagnetic-ferromagnetic-antiferromagnetic 
Ising chain. 

Replacing all Heisenberg operators in Eq. ([]]) by the Ising variables s, t, u one can obtain the 
model, describing by the follow Hamiltonian: 

- f3H = ^(Jf(sA + Uui) - J A UiS i+1 + h(si + U + Ui)). (3) 

i 

Here the Jp = (3Jf and J a = 13 J a are dimensionless ferromagnetic and antiferromagnetic 
constants respectively, h = f3H is dimensionless magnetic field and all the variables s, t, u take 
values equal to ±1. The partition function for one dimensional systems can be written in the 
following way: 

Z = e ~ m = Ylexp(J F (siU + Uui) - J A UiS i+ i + h(si + U + m)). (4) 

(s,t,u) (s,t,u) i 

Here the product is going over all sites of the chain. If the periodic boundary conditions are 
assumed, the partition function can be represented as a trace of product of the transfer- matrices: 

Z = Sp(T N ). (5) 

In our case we can write down the elements of the transfer-matrix as 



T s s = ^2 exp( Jp(st + tu) — Jafus + h(s + t + u)) = ^ exp{a\{t)s + ct2(u)s + as(t, u)), (6) 

t,u t,u 







where the following notations are introduced: 



a\ = J Ft + h, oii = —Jau, «3 = Jput + h(t + u). (7) 

Implementing the summation on t and u in Eq. (^) one can easily obtain the explicit form of 
the transfer-matrix: 



2e jA cosh h + 2e~ jA+2h cosh(2 J F + h) 2e~ jA cosh h + 2e jA+2h cosh(2 J F + h) 
2e~ jA cosh h + 2e jA ~ 2h cosh(2 J F - h) 2e jA cosh h + 2e~ jA ~ 2h cosh(2 J F - h) 



(8) 



Following the general procedure, using the properties of the matrix trace, Sp(Ai A 2 . . . AkB) = 
Sp(BA!A 2 . . . Ak) it is easy to show that the partition function from Eq. (^) for one- 
dimensional model is the sum of the N-th power of the eigenvalues of transfer-matrix. 

z = Y,K N - (9) 

i 

It is easy to see that in the thermodynamical limit, iV — * oo, only the maximal eigenvalue 
contributes to the r. h. s. of Eq. Thus, the problem of finding the partition function 

for the chain reduces to the problem of finding the maximal eigenvalue of the transfer-matrix. 
Solving the secular equation for the matrix (H), one can obtain the eigenvalues 

A± = 2 cosh /i + e" Ja (e 2h cosh(2 J F + h) + e" 2/l cosh(2 J F - h) ± (10) 
± ({e 2h cosh(2 J F + h)- e~ 2h cosh(2 J F - h)) 2 



+ 4(cosh/i + e 2{jA+h) cosh(2J F + h))(coshh + e 2{jA ~ h) cosh(2J F - h))^ ' ' 

It is easy to see that for our purpose we have to take A + . 

Thus, in the thermodynamical limit the free energy per site takes the form: 



N 

N-^o \ 3Nfi 



lim l- lQgA t 1= (111 



~ log 1 2 cosh h + e~ jA (e 2h cosh(2 J F + h) + e' 2h cosh(2 J F - h) + 



3p 

+ ((e 2h cosh(2 J F + h) - e~ 2h cosh(2 J F - h)f 



+ 4{ cosh h + e 2{jA+h) cosh(2 J F + h))(coek h + e 2{jA ~ h) cosh(2 J F - h)}) 1/2 

The factor 1/3 is introduced here because there are three spins in each site of the chain, and 
the total number of spins is 3N. Now, using the general thermodynamical relations we can find 
the magnetization per site 

(dF\ 1 (d\+\ 1 (d\+\ 

The explicit form of this expression is too cumbersome and we found no reason to write it down 
here. However, it is quite suitable to plot analytically the corresponding magnetization profiles 
(diagrams of m versus the external magnetic field H) with the aid of MATHEMATICA 4.0 
package. The results are presented in the next section. 



3 Magnetization Processes. 



It is easy to see that the ground state of the system under consideration at T = and H = 
is the antiferromagnetic spatially modulated structure in which the trimers of spins pointed 
up are alternating with the trimers of spins pointed down (. . . TTtllltttill • • • an d so on). 
Usually this modulated phase is denoted by (3) PTfl . This is valid for all values of J a and Jf 
in the absence of the external magnetic field []. 

The profile of the magnetization as a function of the external magnetic field essentially 
depends on the ratio of coupling constants k = J \j Jp and temperature. Actually one can 
distinguish two regions of values of k: strong antiferromagnetic coupling and weak antiferro- 
magnetic coupling. The limiting case of the latter, when the system reduces to the ordinary 
ferromagnetic Ising chain is J a = 0. When the antiferromagnetic coupling J a is included then 
the low temperature behaviour changes steeply. At T = and H = as is mentioned above 
the chain is in the ordered state with complex antiferromagnetic long-range order. At finite 
but sufficiently low temperatures the long-range order is breaking just by few solitons. The 
external magnetic field will have only small effects until it reaches a critical value H c where it 
can flip a spin without cost in energy. Thus, the appearance of a magnetization plateau in the 
magnetization versus the external magnetic field curve is expected in this case. The width of 
the plateau is obviously equal to H C1 . For the chain at T = one, using the simple principles 
of minimal energy, can easily express H Cl in terms of J a and Jp. 

Let us first consider the case when antifferromagnetic coupling between trimers is assumed 
to be so intensive that under the effect of the sufficiently low external magnetic field only the 
central spins of ferromagneticaly coupled trimers flip (do responde); whereas the rest two spins 
are maintained in the same position as in the ground state. This regime we will call the strong 
antiferromagnetic coupling. So, in this case, when the magnitude of the external magnetic field 
reaches its critical value H Cl the system will pass from its ground state (3) to the novel spatially 
modulated structure (3111), in which the periodic sequence of spins consists of one trimer of 
spins pointed along the field and another trimer with alternating orientation of spins (See Fig. 
1). In order to obtain the value of the critical field H Cl one must compare the energies of (3) 
and (3111) states and find the value of H at which the (3111) state becomes energetically more 
preferable than the ground state. Considering the single periodic spin sequences we have 

- J A - 2H < -AJ F - J A , or H Cl = 2J F . (13) 

So, for the strong antiferromagnetic coupling regime the value of the critical field does not 
depend on J a- 

Now we can give a quantitative criterion for this domain of parameters. The antiferromag- 
netic interaction between edge spins of trimers is strong in the sense described above when the 
flip of the central spin does not lead to the flip of its adjacent one at the critical value of the 
external field calculated above. In the other words, the values of coupling J a and Jp are such 
that the (3111) structure is energetically more preferable than the spin configuration where 
at least one of the spins adjacent to the central one is pointed up at H = H Cl . Comparing 
corresponding energies for the single periodic sequence one can obtain the following condition 
of strong ferromagnetic coupling (Fig. 1): 

1 Hereafter for convenience we will identify J a by J a and Jf by Jf- 
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- 2.1 A - 2H Cl < -2.1 F - 4H C1 , or J A > 3J F . (14) 

The weak antiferromagnetic coupling, Ja < 3Jf, supposes that the ferromagnetic inter- 
action inside the trimers is strong enough, so that the three spins the trimer consists of be- 
have like single Ising spin, which can take values ±3 with effective antiferromagnetic coupling 
Jeff — J a/ 3. This case corresponds to the S = 3/2 singlets in the quantum model considered 
by Hida 0. Thus, for Ja < 3Jp the system under consideration can be regarded as the simple 
antiferromagetic Ising chain with the restrictions given above. The value of the critical field 
magnitude can be easily obtained in this case. Obviously, H Cl = U e ff = \J a, thus, 




if J a > 3 Jf /-. r\ 

if J A < 3J F . 1 ' 

The plots, presented in the Fig. 2 illustrate these results. Note, that the high temperature 
magnetic behaviour for all regions of parameters is rather uniform and the typical magnetization 
profile is given by the curve of Langevin type. So, in that follows we will consider only low 
temperature region. At finite temperatures the values of H Cl are obviously different from Eq. 



(15 ) due to existence of thermal fluctuations of spins. This difference becomes smaller with the 



decrease of temperature. So, for sufficiently low temperatures 

H C1 (T) = H C1 -6(T), hm0(T) = O. (16) 

Fig.2(a)-2(c) show that already at T = 0.1 Jp value of critical field is quite close to \Ja- 

The regime of strong antiferromagnetic coupling is characterized by additional feature which 
is absent in case of k < 3. It is the magnetization plateau at m/m sat = 1/3. Corresponding 
plots are presented in Fig. 3. One can easily see that the width of plateau depends on the 
value of k, namely it becomes broader with the increase of k. In this case another critical value 
of the external magnetic field magnitude defining the width of the plateau appears. So, in the 
magnetization curve one can see a horizontal region for the values of H in interval [H C1 ,H C2 \. 
With the aid of calculations analogous to those we made to obtain H Cl one can easily get the 
second critical value H C2 for the chain at T=0 in the strong antiferromagnetic coupling region. 
Appearance of the plateau at m/m sat = 1/3 and H = H Cl is indication of the stability of (3111) 
configuration when H e [H C1 , H C2 ] . So, when H reaches its second critical value H C2 the above 
mentioned structure begins to destroy due to further flip of the spins which previously were 
pointed down. The comparison of the corresponding energies of single periodic spin sequence 
yields (Fig. 1): 

— 2 Jp — AH < —2.1a — 2H, or H C2 = J A -J F , (17) 

this is in full agreement with Fig. 3. The effects of finiteness of temperature play the same 
role as in the case of H C1 i. e. they can only bring a slight decrease of the H C2 at low values 
of T. It is worth to note once more that high temperatures always lead to disappearance of 
plateau due to strong intensity of thermal fluctuations of spins which in their turn destroy 



the ordered configuration. From all stated above we can conclude that the ferromagnetic- 
ferromagnetic-antiferromagnetic trimerized Ising chain in the region of strong antiferromagnetic 
coupling have a magnetization plateau, its width in the limit of zero temperature is equal to 
H C2 — H Cl = J a — 3 Jf. 

In the Fig. 4 the plots of magnetization for given values of coupling (k = 6) and different 
values of temperature are shown. At T = 2.5 Jf the curve has no deviations from the usual high 
temperature magnetic behaviour of spin systems (Fig. 4(a)). But with decrease of temperature 
some peculiarities appear in the course of magnetization curve. In Fig. 4(b) one can see the 
first signs of forthcoming plateaus at m = and m/m sat = 1/3 and the complete formation 
of these plateaus at T = O.Up (Fig. 4(c)). Fig. 4(c) also clearly indicates that the values 
of critical fields are just slightly different from H Cl and H C2 calculated above. And finally the 
curve, presented in Fig. 4(d), obtained for extremely low temperature T = 0. 001 Jf shows 
step-like behaviour which is typical for T = 0. In this case, one can see that the jumps in the 
course of the magnetization exactly coincide with H = Up and H = J a — Jf that is in full 
agreement with the pattern of transition from (3) state to (3111) state with further transition 
to the saturated state. Nevertheless, it is necessary to emphasize that there are no real jumps 
at any finite temperatures at all, because the system under consideration is one-dimensional 
and it is know that there are no phase transition in classical systems in d = 1. So, the results 
presented in Fig. 4(d) should be regarded not as jumps of magnetization but as extremely 
steep rises of the m(H) function in the vicinity of critical values of H. 

4 Conclusion. 

The main aim being pursued in this paper is to establish the formation of magnetization 
plateaus in the classical one-dimensional spin systems with Ising interaction. As an example 
we chose the classical counterpart of the ferromagnetic-ferromagnetic-antiferromagnetic trimer- 
ized Heisenberg chain, considered by Hida . The approach developed by Hida was based on a 
numerical diagonalization of finite size system and allowed him to explain the low temperature 
magnetization data for the compound 3CuCl2 ■ 2dx as well as to predict the appearance of a 
magnetization plateau at m/m sa t = 1/3 when the ferromagnetic exchange energy is comparable 
to or smaller than the antiferromagnetic exchange energy. Using entirely analytic technique of 
the transfer matrix for corresponding system of Ising spins we have obtained the exact explicit 
expression for magnetization function. Despite the extremal complexity of the expression, nev- 
ertheless it permits to draw the plots of magnetization processes for arbitrary finite temperature 
and couplings. Obtained plots completely reproduce the result of Hida, moreover in our model 
the exact condition of existence of plateau J a > 3Jf is derived. It is also possible using our 
approach to calculate the width of the plateau at m/m sat = 1/3. The latter is equal to J4 — 3Jf 
in the limit T — > 0. 

The phenomenon of magnetization plateaus are often regarded to have purely quantum 
origin, however it is obviously not always true since there are some two dimensional classical 
spin models which exhibit a magnetization plateaus at zero temperature. For instance, in Ref. 
[|T2"fl Kubo and Momoi investigated the Multiple Spin Exchange (MSE) model on the triangular 
lattice in the classical limit and have found a large range of parameters where magnetization 
plateaus at m/m sat = 1/3 and m/m sat = 1/2 occur at zero temperature. There are also results 



on the two dimensional Ising models [16| as well as on the calculations on recurrent lattice [20 
However, for one dimensional systems these results are obtained for the first time. 

Moreover, these allow to suppose that there exists a class of one dimensional spin systems 
which being both considered in terms of Heisenberg and Ising spins leads to the qualitatively 
same structures of the magnetization profiles, particularly to the formation of magnetization 
plateaus. In contrast to the quantum systems the plateau in our model has very simple nature 
associated with the stability of spatially modulated structure, so called (3111) configuration of 
spins which we proved to exists at T = 0. Apparently, for finite but sufficiently low temperatures 
the plateau is still caused by the same (3111) structure, though the other configurations are 
also possible, for instance so called uud phase which is responsible for plateau at m/m sat = 1/3 
in triangular antiferromagnets and MSE model. 

In conclusion note, that the approach based on Ising spins is remarkable for its simplicity 
and clarity. It doesn't require any numerical calculations in contrast to the majority of methods, 
which are used in this field of research, and can serve as a very useful supplement to the generally 
applicable techniques such as the numerical diagonalizations of finite clusters, Quantum Monte- 
Carlo calculations, series expansions and bosonization technique. 
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Figure Capture 



Fig. 1 Periodic spin sequences and corresponding energies for (a): ground state configu- 
ration at T = and H = 0; (b): the plateau state, so called (3111) configuration; (c): one of 
the possible intermediate configurations which forms after the plateau state. 
Double line between the sites denote the antiferromagnetic interaction J a- 

Fig. 2 Magnetization curves in the regime of weak antiferromagnetic coupling, Ja < 3Jp 
at T = O.lJp for different values of k: (a) k — 1; (b) k — 2; (c) k — 3. 
The value of the critical field magnitude in these cases is 3/2 J a- 

Fig. 3 Magnetization curves in the regime of strong antiferromagnetic coupling exhibiting 
the magnetization plateau at m/m sat = 1/3 for k — 4 (a), k — 5 (b) and k — 10 at T — O.lJp. 
The values of critical field magnitude, determining the width of plateau are H C1 = 2J F and 
H C2 = Ja — Jf- 

Fig. 4 The magnetization curves for k = 6 at different temperatures: (a) T = 2.5J/?; (b) 
T = 0.8J F ; (c) T = 0.1 J F ; (d) T = 0.001J F . 
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